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Abstract. We predict the charm and bottom quark nuclear modification factors 
using weakly coupled pQCD and strongly coupled AdS/CFT drag methods. The 
\og{pT / Mq) / pt dependence of pQCD loss and the momentum independence of 
drag loss lead to different momentum dependencies for the Raa predictions. This 
difference is enhanced by examining a new experimental observable, the double 
ratio of charm to bottom nuclear modification factors, R^^ = R-aa/ ^AA- LHC 
the weakly coupled theory predicts R^^ 1 whereas the strongly coupled theory 
predicts R*^^ ~ .2 independent of pT- At RHIC the differences are less dramatic, 
as the production spectra are harder, but the drag formula is applicable to higher 
momenta, due to the lower temperature. 
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1. Introduction 

The excellent agreement between the pion nuclear modification factor, R\j^{pt)^ 
and the perturbative quantum chromodynamics (pQCD) prediction, combined with 
the null control experiments of d + A, direct photons, and the consistency check of 
R\j^{Pt) ^ R\a(Pt) [1., makes a compelling argument for the creation of a weakly 
coupled quark-gluon plasma (QGP) at RHIC. However this picture of a weakly coupled 
plasma is undermined by the failure of pQCD to simultaneously account for several 
other experimental observables: the nearly ideal hydrodynamic fiow of low-p^ particles 
(requiring an entropy to viscosity ratio rj/s ~ .1); the highly suppressed (similar 
in magnitude to the pions) nonphotonic electrons, the decay fragments of heavy 
charm and bottom quarks; and the nonperturbatively large intermediate-pT azimuthal 
anisotropy, V2 [2 . Much progress has been made in reducing the discrepancies between 
pQCD and data; one cannot currently claim the falsification of pQCD techniques 
applied to RHIC Nonetheless, one might naturally consider the possibility that 
strongly coupled dynamics might provide a better theoretical description of data. 

Exploiting the AdS/CFT correspondence, which allows many strongly coupled 
supersymmetric field theoretic problems to be solved analytically in a classical 
supergravity (SUGRA) dual, a number of calculations have been made that are, in 
fact, in qualitative agreement with known results. The infinitely coupled entropy 
density is three-fourths the value of the Stefan-Boltzmann limit, similar to that seen 
on the lattice. There appears to be a universal lower bound of entropy to viscosity. 
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v/^ ^ l/47r ~ .1. Mach wave-like structures appear for supersonic motion. And heavy 
quarks moving through strongly coupled plasma appear to lose a significant fraction 
of their momentum [4 . Despite these successes, super Yang-Mills (SYM) is quite 
different from QCD and, partly for this reason, strong coupling calculations have so 
far failed to make precision predictions. 

What we seek, then, is a robust, distinguishing measurement that could possibly 
falsify one or both approaches. It turns out that measuring the heavy quark energy 
loss for charm and bottom quarks separately provides just such a probe. 



2. AdS/CFT Energy Loss Compared to pQCD 

The AdS/CFT conjecture has been applied to three types of energy loss models. The 
first uses the Wilson line formulation for the radiative transport coefficient q; the 
conjecture allows a numerical estimate of its value. The second finds the momentum 
diffusion coefficient, used in a collisional heavy quark energy loss computation. 
The third models the entire plasma plus heavy quark system in the AdS/CFT dual 
and derives the drag coefficient [5] . There is considerable ongoing debate and research 
on these approaches (see, e.g. [5l, ^61); we will focus on the latter because it does not 
use the AdS / CFT correspondence to find the value of an input parameter needed in a 
pQCD-inspired model. Rather its more ambitious goal of modeling the entire energy 
loss process in the SUGRA dual leads to novel mass and momentum dependencies. 

The drag on a heavy quark of mass Mq moving with constant velocity in a 
constant temperature SYM plasma in the limit of A = \/ 9sym^c ^ 1, A/'c ^ 1, 
Mq > T^^^ is [5^ 

^ = -^^^^ = 2Mg 

Issues related to the relaxation of the strong assumptions made in deriving Eq. ([T]) and 
the momenta above which one begins to expect significant corrections are discussed 
later in the text. 

One typically does not consider drag, but rather fractional energy loss e, where 
Pf = (1 — e)pi. For approximately constant, luminal velocity the average fractional 
loss for a quark propagating through a medium of length L is CAdS = 1 — exp(— /igL), 
independent of pt- On the other hand the asymptotic form of the average pQCD 
radiative fractional energy loss is CpQCD = '^(^sQL'^ log{pT /Mq) /pt, where is a. 
dimensionless proportionality constant. In this case I^qcd ^ as increases. 

The exceptional momentum reach of LHC and its particle identification (PID) 
capabilities can be used to distinguish between these two. The production spectrum 
for heavy quarks is well approximated by a power law, dNQ/dpT - 1/pt'^^^^\ with 

driQ/dpT > 0. We then have R^AiPT) = ((1 — ^))^^*^^^^~"^)^eom, where, in the full 
numerical results shown in the figures the average over realistic geometry includes jet 
production proportional to the binary distribution and propagation through a medium 
whose transverse profile is proportional to the participant distribution. Bjorken 
expansion is included, and a diffuse Woods-Saxon base nuclear profile is used. nQ{pT) 
was found from a best fit of FONLL spectra [7 . 

For LHC the increase in the power law index is not fast enough to compensate for 
the decrease in pQCD fractional energy loss; thus dBF^^^ /dpT > 0. Contrariwise, the 
increase in the index combined with the momentum independence of the AdS/CFT 
drag implies dR^j^ /dpr < 0. 
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To compare the AdS/CFT drag predictions to pQCD predictions and future 
experimental data one must connect SYM to QCD. This requires a nontrivial mapping 
of the running couphng in QCD to the constant couphng in SYM and of the 
temperatures. The most naive we denote as "obvious"; this takes asvM = <^s = const. 
and T^^^ = T^^^ . An "alternate" might equate the quark- ant iquark force from 
AdS/CFT to that measured on the lattice, yielding Xsym — 5.5, and the energy 
densities, yielding, due to the approximately factor of three greater number of degrees 
of freedom in SYM, T^^^ = T^cd ^ 

As an illustration of one of the difficulties in making a precision strong coupling 
prediction consider the momentum diffusion coefficient mentioned earlier, D = 
2l\f\KT 0. Using the "obvious" prescription and ots = .3, a common value taken 
in pQCD jet physics [3], 1.2/27rT. However D :^ 3/27rT was claimed to better 

reproduce experimental data [5]. This would require a coupling of ag .05, far from 
the strong coupling limit. Nevertheless we will use these two a^'s as representative 
values when employing the "obvious" prescription. We note that a smaller coupling 
leads to a larger diffusion coefficient because D is in the coordinate representation 
instead of the usual momentum one. 

Eq. ([T]) is implemented numerically as in \J_^. Energy loss is assumed to start 
at thermalization, tq ~ .6 — 1 fm/c, and stops when the confinement temperature, 
Tc ~ 160 MeV, is reached. The exponentiated dependence in /ig leads to an 
extremely sensitive probe of the opacity of the medium, as well as tq and Tc. We will 
employ the WHDG model of convolved elastic and radiative energy loss for pQCD 
predictions |3j . LHC predictions using the ASW model of purely radiative energy loss 
have also been made with extrapolations of extreme 40 < g < 100 [8]; instead of using 
the exact ASW model we approximate these predictions with radiative only WHDG. 
In all calculations initial state effects were neglected; just as at RHIC for definitive 
experimental statements to be made a p + Pb control run will be essential. 
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Figure 1. (Color Online) Raa{pt) and R\a(jPt) for central Pb + Pb 
reactions at LHC with energy loss given by AdS/CFT drag and by pQCD 
using WHDG [3 . The full numerics produce the expected behavior of 
dR^j^ /dpT < and dR^^^^ /dpr > 0. Representative possible initial 
gluon densities are given by dNg/dy = 1750 and dNg/dy = 2900 from 
the PHOBOS extrapolation and the KLN model of the CGC, respectively 
[9]. The AdS/CFT drag parameters are from the "obvious" prescription 
with a SYM = .05, giving D = 3/27rT (abbreviated to D = 3 in the figure). 
The "O" and "|" symbols come from Eq. Q and are explained in the text. 

In Fig.[l] we see that our expectations for the momentum dependence of the charm 
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and bottom Raa for the two theories are borne out in the fuh numerical calculations. 
While not shown, one is not surprised that the purely radiative, highly suppressed, 
large q extrapolations have momentum dependencies difficult to distinguish from some 
of the AdS/CFT drag predictions. It turns out that we can enhance the signal by 
taking the double ratio R^^{pt) = R^a^Pt) / R\a{Pt) - From the asymptotic fractional 
energy loss formula above pQCD gives 



Kqcd{Pt) ^ 1 - ^asn{pT)L' \og{Mi,/M,) 



Pt 



(2) 



where ni){pT) — ndpr) — ^{pt)^ and the argument of the log is the ratio of bottom to 
charm quark mass. Hence pQCD predicts that this ratio goes to unity asymptotically, 
and that this approach is slower for greater suppression. 

Taking the medium to be have a uniform, time-independent density of size L, 
one may estimate the AdS/CFT drag nuclear modification factor: Raa ^ ^/^Ql^Q^- 
Then 



RAdsM 



rihlJih 



(3) 



riclJic Mb 

and AdS/CFT drag predicts R'^^ ~ .27, independent of momentum. 

One sees in Fig. [2] that the above approximations are well reproduced by the 
full numerical results. Most normalization differences for the wide range of input 
parameters used cancel, and the curves bunch into a pQCD group and an AdS/CFT 
drag group. Additionally the more highly quenched pQCD curves approach 1 more 
slowly, except for the highly oversuppressed q = 100. Supposing that LHC data are 
similar to the pQCD predictions one might be able to distinguish between convolved 
elastic and inelastic loss, for which R*^^ monotonically increases, and purely radiative 
energy loss, for which R^^ dips to a minimum near pt ~ 10 GeV; one must be cautious, 
though, as many of the assumptions in the elastic energy loss derivations break down 
for Pt ^ Mq. 
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Figure 2. The double ratio of Raa{Pt) to R!aa{Pt) predictions for LHC 
using Eq. ^ for AdS/CFT and WHDG [3] for pQCD with a wide range of 
input parameters. By taking the ratio the large normalization differences 
cancel, and the results collect into two distinct groups. 

It turns out, though, that there is reason to believe that Eq. ([T]) cannot be applied 
to arbitrarily large pt- For strictly finite heavy quark mass, demanding a time-like 
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string endpoint in constant velocity motion gives a "speed limit" of [TO] 



2M 



4M2 



One may arrive at this same value through a different line of reasoning. Assume that 
the quark's constant velocity is maintained by an electric field. The largest electric 
field sustainable from the Born-Infeld action limits the magnitude of the momentum 
loss; the critical speed after which the field cannot be strong enough to keep the quark 
velocity constant is given by Eq. [10 . 

Eq. Q was found assuming a constant plasma temperature; this is not the case 
in experiment. To give some idea of the momenta above which one might expect 
significant corrections to Eq. ([T]) we plot the smallest possible 7c (corresponding 
to the largest temperature, T{x = 0,t = tq) and the largest 7c (from the smallest 
temperature, Tc). These are indicated by an "O" and a "|" in the figures, respectively. 

Future RHIC detector upgrades will allow for individual charm and bottom quark 
detection. Predictions for R\j^{pt) and R\j^{pt) at RHIC from Eq. ([l]) and pQCD 
are shown in Fig. |3] (b). As seen in Fig. |3] (a) the power law production index grows 
quickly at RHIC; we no longer can expect that the decrease in pQCD energy loss as a 
function of momentum will result in dEF^^^ /dpT > 0. In fact Fig. [s] (b) shows that 
the full numerical results for from pQCD and AdS/CFT drag both decrease with 
PT' Nonetheless one may still examine the double ratio i?^^. Fig. [ij While the larger 
index makes the grouping less dramatic at RHIC one may still differentiate between 
pQCD and AdS/CFT drag. Due to its smaller multiplicities, the temperature of the 
medium at RHIC is smaller than will be seen at LHC; hence the AdS/CFT drag 
"speed limit", Eq. (|4|, is higher at RHIC than LHC. 




— RHIC charm 
■- RHIC bottom 

— - LHC charm 

— ■ LHC bottom 



0.8 



3 



Charm - 

I 

1 pQCD Rad+EI 

J -- AdS/CFT, D= 3 
■-■ AdS/CFT, D=1 
-j --- AdS/CFT, ^ = 5.5 

1} ; 


1 Bottom 

- A 







50 100 

(a) 



20 30 
P^(GeV) 



50 
RHIC 

(b) 



20 30 

p,(GeV) 



Figure 3. (a) The power law production index nQ{pT) for RHIC and 
LHC. While it is quite flat over a large momentum range at LHC, at RHIC 
^q(pt) hardens appreciably as momentum increases, (b) Raa{pt) and 
Raa{pt) for RHIC using AdS/CFT drag and pQCD WHDG. The large 
increase in uq^pt) overcomes the decreasing fractional momentum loss for 
pQCD; both AdS/CFT drag and pQCD results decrease as a function of 
momentum at RHIC. 



3. Conclusions 



We compared the R^j^{pt) predictions from weakly coupled pQCD and strongly 
coupled AdS/CFT. The momentum dependencies of the energy loss formulae resulted 
in dR^^^^ /dpT > and dR^^ / dpr < at LHC. The difference in the mass and 
momentum scales for the two make the new observable R^^{pt) = R\j^{pt) / R\a{Pt) 
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Figure 4. The double ratio of Raa{pt) to R\a{pt) predictions for RHIC 
using Eq. ^ for AdS/CFT and WHDG [3 for pQCD with a range of input 
parameters. While the hardening of the production spectrum reduces the 
dramatic bunching at RHIC as compared to LHC, the lower temperature 
at RHIC means the AdS/CFT drag formalism is applicable to higher 
momenta. Note that R""^ is plotted to only 50 GeV for RHIC. 

an excellent probe for falsifying a coupling limit. These calculations were also 
performed for RHIC. The harder power law index meant that both pQCD and 
AdS/CFT drag predicted decreasing dR^^/dpT < 0. The double ratio R^^, though, 
still had differentiating power. Two issues are of concern: momentum range of 
applicability and universality. pQCD calculations may only be valid at very high 
Pt] AdS/CFT drag calculations may only be valid up to some lower momentum limit. 
It is not impossible that an overlap region for which both ought to be valid does 
not exist. The possibility of falsification rests with unique, measurable differences in 
theoretical predictions. A theory which does not give this may not be falsified by the 

data, but is also of little use. 
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